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We discuss nontrivial features of the large scale structure of the universe in the simplest curvaton 
model proposed in our paper astro-ph/9610219 The amplitude of metric perturbations in this model 
takes different values in different parts of the universe. The spatial distribution of the amplitude 
looks like a web consisting of exponentially large cells. Depending on the relation between the cell 
size Ao and the scale of the horizon Ih, one may either live in a part of the universe dominated by 
gaussian perturbations (inside a cell with Ao S> Ih), or in the universe dominated by nongaussian 
perturbations (for Ao ^ Ih)- We show that the curvaton contribution to the total amplitude of 
adiabatic density perturbations can be strongly suppressed if the energy density of the universe 
prior to the curvaton decay was dominated not by the classical curvaton field but by the curvaton 
particles produced during reheating. We describe the curvaton-inflaton transmutation effect: The 
same field in difi'erent parts of the universe may play either the role of the curvaton or the role of the 
infiaton. Finally, we discuss an interplay between the curvaton web and anthropic considerations in 
the string theory landscape. 
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I. INTRODUCTION 



It is well known that quantum fluctuations of the in- 
fiaton field in the simplest inflationary models produce 
gaussian adiabatic perturbations of metric with nearly 
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scale independent spectrum 0, 0, ■ This prediction is 
in a verygood agreement with the present observational 
data 0, U- However, since the precision of the observa- 
tional data improves every year, one may try to discrim- 
inate between various inflationary models by considering 
subtle deviations from the standard paradigm. One way 
to distinguish between various models is by studying a 
small degree of nongaussianity of perturbations, and a 
possible small contribution of the isocurvature perturba- 
tions. 

It is impossible to produce isocurvature perturbations 
in the models describing a single scalar fleld, the inflaton, 
and the degree of nongaussianity of the inflaton pertur- 
bations typically is extremely small. However, it is rela- 
tively easy to produce isocurvature and/or nongaussian 
adiabatic perturbations in the models containing several 
light scalar fields. The first model of that type was pro- 
posed in It was based on investigation of axions, 
which were massless during inflation. At that stage, the 
axion perturbations practically did not contribute to the 
perturbations of curvature. At the late stages of the evo- 
lution of the universe, the axion field became massive, 
and its initial isocurvature perturbations were converted 
into a mixture of adiabatic and isocurvature perturba- 
tions, which could also have a significant nongaussian 
component. Another model of that type describing a 
mixture of adiabatic and isocurvature perturbations was 
considered in 0. 

However, current observational data put very strong 
constraints on thepossible admixture of isocurvature per- 
turbations 0, 0, 111 ■ Therefore recently a special atten- 
tion was payed to a class of models where a second field, 
called "curvaton," at a certain stage begins to dominate 
and then decays leading essentially to a secondary stage 
of reheating. In this model the initial isocurvature fluc- 
tuations can be completely converted to the adiabatic 
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ones. The spectrum of perturbations in such models can 
be nonfiat, it can have various nonstandard features, and 
it can be nongaussan. This possibihty was first discussed 
in our paper [9|, and then it was significantly developed 
in the series of papers by Lyth, Wands lOj, Moroi and 
Takahashi fn] and other authors. A closely related but 
different mechanism was also proposed in 

At the first glance, our original model may seem ex- 
tremely simple. One of its versions describes a light field 
(T, which fluctuates during inflation, and later begins to 
dominate the energy density of the universe and decays. 
However, the simplicity of this model, as well as of many 
subsequent versions of the curvaton scenario, is some- 
what deceptive. Such models must satisfy many speciflc 
conditions on the mass of the curvaton, on the way it de- 
cays, on the way the CMB and the baryon asymmetry of 
the universe are produced, etc., see e.g. iS^ Ji)j JjJ. Also, 
the predictions of the curvaton scenario may depend on 
the choice of initial conditions and on our own position 
in the universe. This last fact, which was mentioned in 
0, will be the main subject of our paper. 

Most of the papers on the curvaton scenario discuss 
the behavior of the curvaton perturbations as a func- 
tion of the classical nearly homogeneous curvaton field. 
Meanwhile we would like to concentrate on the global 
structure of the universe in this scenario, in which the 
curvaton field may take different values in different pars 
of the universe due to stochastic processes during infia- 
tion. We will show that the spatial distribution of the 
amplitude of density perturbations in the curvaton sce- 
nario looks like a web consisting of exponentially large 
cells. We will call this cell structure 'the curvaton web.' 
Depending on the values of the parameters of the theory, 
the amplitude of density perturbations reaches its maxi- 
mal values either in the internal parts of the cells, or near 
their boundaries. 

We will also describe a potentially important effect, 
which, to the best of our knowledge, was ignored in the 
existing literature on the curvaton perturbations. It is of- 
ten assumed that the curvaton energy density is entirely 
due to the initial classical curvaton field and its infia- 
tionary perturbations, which also behave essentially as 
classical fields. However, part of the infiaton energy can 
be converted to the energy of curvaton particles during 
reheating after infiation. We will show that under cer- 
tain conditions the energy density of these particles can 
be much greater than the energy density of the primor- 
dial curvaton field. As a result, this effect may reduce the 
curvaton contribution to the total amplitude of density 
perturbations. It may also allow to increase the degree 
of nongaussianity of these perturbations. 

Finally, we will discuss the curvaton-infiaton transmu- 
tations and a possible interplay between the curvaton 
web and anthropic arguments related to the concept of 
string theory landscape. 



II. GENERATION OF PERTURBATIONS 

Consider a theory of a scalar field a with a mass to, 
which gives a subdominant contribution to the energy 
density during infiation. If m is much smaller than the 
Hubble constant H during infiation, the long-wavelength 
perturbations of this field with an average amplitude 
5(J = are generated during each time interval H~^. 
Ignoring the spatial variations of the classical field cr, 
which are locally unimportant but can be significant on 
a sufficiently large scale (see below) , one finds that these 
perturbations give rise to locally gaussian perturbations 
of density 5pa — m^aSa ~ "'^^'^ . Their relative ampli- 
tude compared to the (local) value of the energy density 
of the field a is given by 

Sp _ 25a H 

Pa (J TTCT 

This simple equation is valid only for da > a; ^ become 
0(1) for a < H/2tt. 

The curvaton scenario is based on the assumption that 
the perturbations — can be significant. However, since 
this field is subdominant during infiation and shortly af- 
ter it, the contribution of the curvaton ffuctuations to 
the adiabatic perturbations of metric related to is 

Ptotal 

much smaller than — : 

dp m^aSa m?a 

Ptotal Ptotal 67riJ 

where we use the system of units = (SttG) ^ = 1. 

After infiation, the infiaton field decays and produces 
matter. If the products of its decay are ultrarelativis- 
tic particles, then their energy density subsequently de- 
creases as a~^, where a is the scale factor of the universe. 
Meanwhile, if the field a has a quadratic potential, then 
for a long time its energy density decreases more slowly. 
It does not decrease at all until the Hubble constant be- 
comes smaller than to, and then it decreases as a~^. As 
a result, the field a may give the dominant contribution 
to the energy density of the universe, and its perturba- 
tions may give rise to significant perturbations of metric. 
If cold dark matter and baryon asymmetry are produced 
during a secondary stage of reheating due to the curva- 
ton decay, then the adiabatic perturbations of density are 
generated, without any admixture of isocurvature pertur- 
bations 1^, Il3 . These perturbations are of the order . 
In principle, they can be much greater than the adiabatic 
perturbations existing at the end of infiation. 

Two comments are in order here. First of all, parti- 
cles a can be produced during the infiaton decay at the 
first stage of reheating. After the moment when they be- 
come nonrelativistic, their energy density decreases in the 
same way as the energy density of an oscillating homo- 
geneous field (T, and of the long- wavelength fiuctuations 
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of the curvaton field. This is an important and rather 
generic effect, which, to the best of our knowledge, was 
not discussed in the literature. It can increase the to- 
tal energy density of the curvaton field at the time of its 
decay, and decrease the resulting amplitude of the adia- 
batic perturbations produced by the decaying curvaton 
field, see Section Hvl 



Another comment concerns the origin of the classical 
scalar field a in the denominator of Eq. This field 

can be nearly homogeneous in the whole universe, with 
the value determined by initial conditions at the begin- 
ning of inflation. However, if the mass of the field is much 
smaller than the Hubble constant during inflation, more 
interesting situation can be realized. Namely, the fluc- 
tuations of the curvaton field produced during (eternal) 
infiation may bring the light field a to different places in 
different parts of the universe. As a result, the quantity 
^ and the resulting amplitude of adiabatic perturba- 
tions of metric take different values in different exponen- 
tially large parts of the universe. This effect, as well 
as the similar effect considered in [14] in the context of 
the inflationary Brans-Dicke cosmology, provides a possi- 
ble justification to the anthropic considerations involving 
variable amplitude of perturbations of metric TB^ . 

In this paper we will describe a simplest version of the 
curvaton scenario proposed in (3, which is sufficient to 
explain our main idea. In the beginning, we will assume, 
for simplicity, that inflation occurs for a very long time, 
and during this time the Hubble constant H and the 
effective mass of the curvaton field practically did not 
change. This is a reasonable assumption in new infiation 
and in some versions of hybrid inflation. In the subse- 
quent sections we will consider possible generalizations 
of these assumptions. 

Under the specified conditions, the scalar field cr during 
the long inflationary evolution will gradually roll down 
towards cr = 0, and then it will jump in all possible di- 
rections due to inflationary perturbations. In this respect 
our simplest model 9] differs from the scenario where the 
universe is filled by a homogeneous curvaton field. As we 
will see in Section IVTl a generalized version of our sce- 
nario in fact describes a collection of exponentially large 
domains with all possible initial conditions for the field 
CT. However, if one calculates the average value of the 
field cr over the whole universe in our scenario, one finds 
(cr) = 0. Meanwhile, for (cr^) in the simplest model with 
H = const one finds the well-known Bunch-Davies ex- 
pression 0,111,113: 



Hi 

47r2 



H 



dk 



(3) 



The estimate of the energy density of the fluctuations of 
the fleld a can be obtained as follows: 



167r2 



(4) 



To get an intuitive understanding of this result, one 
should take into account that the fluctuations of the field 
cr with mass ^ produced during each time in- 
terval At = have an average amplitude = ^ 
0, 0, ^3 ■ I'^ beginning, an average square of the 



amplitude of these fluctuations grows as (a'^) 



Ml 
How 



just like in the theory of a massless scalar field 
ever, the amplitude of long- wavelength fluctuations of a 

massive scalar field decreases as a{t) — aoexp (^~Y[r\ 
Therefore the amplitude of the fluctuations will start de- 
creasing exponentially after the time At ^ H/m^ since 
the moment they were produced. As a result, the main 
contribution to (cr^) will be given by perturbations pro- 
duced during the latest time interval At ~ because 
all other perturbations will become exponentially small. 



This gives the estimate (cr^ 



^ = T^l^, which 



almost coincides with the Bunch-Davies result. 

Note that inflationary fluctuations with k <^ H cor- 
respond to quantum fluctuations with large occupation 
numbers Uk — ^ ^ 1 Therefore these fluctuations 
can be interpreted as waves of a classical scalar field. The 
quantity 
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2 27rm 



(5) 



represents an average deviation of the inhomogeneous 
classical scalar fleld from zero on an exponentially large 
scale ^ ^ Ao, to be discussed in the next section. 



III. THE CURVATON WEB 



If we decide, as it is commonly done, to concentrate our 
attention only on (cr) = and (cr^) — ^Ji^i , or on any 
other averages involving an arbitrary number of opera- 
tors cr at different points, we will be unable to get a full 
picture of the distribution of perturbations within rele- 
vant for observations regions. In order to go beyond the 
simple perturbative approach, one should first remember 
that as a result of accumulation of the long-wavelength 
quantum fluctuations of the field cr, the universe becomes 
divided into many exponentially large regions filled with 
the positive field a with the typical value cr ^ ^Yirmi 
and many exponentially large regions containing nega- 
tive field cr with the typical value cr ^ — tt^ — . 

In order to understand the properties of the "curva- 
ton landscape" one should take into account that for 
the observers interested in the perturbations on the scale 
I ~ k~^, all perturbations with the wavelength A ^ k~^ 
will look like a homogeneous field. For these perturba- 
tions, which look locally homogeneous on a scale fc~^, we 
have 



(6) 
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One can easily check that 50% of the total Banch- 
Davies value ^^^i is given by perturbations which look 
relatively homogeneous on scale smaller than A0.5 ^ 

ff^^exp^^^ and 90% of the Banch-Davies value is 
given by perturbations which look homogeneous on scale 



smaller than Xq q '-^ H exp(^^^ j . Note that the scales 

exp^^j and H^^ exp^^^ j are exponentially dif- 
ferent from each other. Therefore it would be incorrect to 
think that the universe after inflation looks like a chess- 
board with cells of positive and negative cr of approxi- 
mately equal size I ~ A0.5 (or I ^ A0.9). Still the loga- 
rithms of the scales A0.5 and A0.9 are similar and propor- 
tional to /m? . In this rather vague sense we will be 
talking about the distribution of the curvaton field being 
locally homogeneous on scale 



Ao 



iJ^^exp 



(7) 



To get a better understanding of the spatial distribution 
of the curvaton field we performed computer simulations 
of the process of formation and accumulation of the cur- 
vaton perturbations, using the methods developed in Ref. 
[13, see Fig.H] 




FIG. 1: Results of a computer simulation of the Bunch- 
Davies distribution of the curvaton field a for H/m — 3. 
For greater values of H/m the typical height of the hills 
grow as /m, and a typical size of the 'islands' grows as 

_H"-'exp( O 



In order to make the level a = clearly 

visible, we show only the part of the 'curvaton landscape' with 
o" > 0, i.e. above the 'sea level' a = 0. As we will show below, 
the amplitude of the curvaton density perturbations behaves 
in a rather peculiar way near the 'shoreline' a = 0, see Figs. 
Hand 13 

Exponentially large domains with positive and nega- 
tive a are separated by the boundaries with a = 0, which 



correspond to the 'shoreline' in Fig. ^ The 3D distribu- 
tion of the domain boundaries with cr = 0, which we will 
call 'the curvaton web,' resembles the network of domain 
walls in the theories with spontaneous breaking of a dis- 
crete symmetry. The main difference is that usually the 
domain walls have large energy density; they exist only 
because of the topological reasons. In our case the situa- 
tion is opposite: The potential energy density of the field 
a is of the order H'^ inside the domains, and it vanishes 
at the boundaries with cr = 0. On the other hand, the 
amplitude of the curvaton density perturbations behaves 
in a rather peculiar way near the walls with cr = 0, see 
Figs. El and 13 

In order to evaluate the perturbations of metric, let us 
first consider perturbations on the scale much larger than 
the size of a typical domain, / ~ ^ Aq 0. In this 
case the fluctuating field cr wanders many times in the 
region —H^/m < cr < H^/m on the scale I, so its value 
averaged over the domain of a size vanishes. As 

a result, addition of a perturbation Sa{k) with k <^ kg 
does not lead to the usual density perturbations w?a6(j 
pl|. Perturbations of density will be nongaussian. To 
get a rough estimate of their amplitude, one may assume 
that they are quadratic in i5cr, which would lead to an 
estimate for 7^ at the stage of inflation. 



m 



(8) 



see Eq. (5) in 9]. (Note that here we are writing 
{pa) = 0{H*) instead of the local value pa{x).) How- 
ever, it was pointed out in our paper j2| that for a more 
accurate evaluation of the amplitude of the perturbations 
one should flnd the correlation function 



1 



{Pa 



Pa{x.t)p„{x + r ,t)) - {p„) 



(9) 



After that one can find 

This yields 
5p 



from the relation 



dk sin kr ( 5p 



k kr 



m^k^ 



{Pa 



k' 



{Po 



fc-fe' 



(Tk 



(10) 



where cr^fc^ = 3 



for k ^ H . Because we 
are only interested in large scale fluctuations produced at 
inflation, we can make a UV cut off at A; ~ H . The corre- 
sponding calculation have been made in Q taking into ac- 
count the complicated inflationary and post-inflationary 
evolution in the chaotic inflation scenario for all wave- 
lengths. However, in the context of the present discus- 
sion we are assuming that the Hubble constant during 
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inflation remains constant, and we will be interested only 
in the long wavelengths which do not loose their gradi- 
ent energy due to the subsequent redshift. In this case 
the calculation is much simpler than the calculation per- 
formed in our paper ^ , but the result is qualitatively the 
same as Eq. (17) of Ref. 0: 



m / k 



(11) 



Note that and the perturbations cr^ scale in the same 
way during the subsequent expansion of the universe 
(they decrease during the field oscillations as a~^), this 
equation remains valid until the perturbations enter the 
horizon. 

The reason why this result is greater than the naive 
estimate (|SJ by a factor of H/m is pretty simple. The 
estimate ^ was obtained by integration of the square of 
perturbations Sa ~ H/2tt. However, the positive (nega- 
tive) perturbations 5a are correlated with the exponen- 
tially large domains where the field a is positive (neg- 
ative) and takes average value ~ /ra, which is H/m 
times greater than 5 a. 

Eq. (|ll|l implies that these perturbations have blue 
spectrum. Note, however, that one can easily obtain red 
spectrum as well, or a spectrum of a more complicated 
shape, if one takes into account a slow decrease of H and 
a possible change of m during inflation. 

We should note also that Eq. (|ll|l describes the ratio 
-gf-^, where {p^) is the average energy density of the field 
(J. On scale I ^ Aq the local energy density p„ changes 
substantially, so in some places, e.g. near the walls of the 
curvaton web, where (7 = 0, the quantities 7^ and — 

\Pcr/ Per 

may be quite different from each other. 

Now let us consider an opposite limit: / <C Aq ^ 

H^^ exp^O^^^^ . On this scale the value of the field a 

remains practically constant, and the relative perturba- 
tions of density of the field a are (locally) gaussian, with 
an amplitude given by Eq. 



5p 



25(7 
a 



H 

TTCT 



(12) 



The last term appears due to the red tilt of the fluc- 
tuations, but this term is indistinguishable from 1 for 
fc-i ~ / < Aq. 

This equation shows that ^ can be large at small a. 

Pa 

Strictly speaking, this equation is not correct in the re- 
gions near cr = 0. Indeed, the field a deviates from ct = 
by (5(7 ± tP- on the scale , which means that the rel- 
ative amplitude of density perturbations becomes large 
near (7 = 0: ^ = (9(1). 

Per 

A typical amplitude of perturbations can be obtained 
by substitution of the average amplitude of a from Eq. 



This gives 

2 27rm 

5p m 



(13) 



Thus, up to a factor (-g-)"^, which is approximately 
equal to 1 for Z <C Aq, the typical amplitude of the fluc- 
tuations on scale I <C Aq (|13|) . even though the nature 
and interpretation of these perturbations is quite differ- 
ent. Meanwhile on scale I 3> Aq the amplitude of pertur- 

bations I|1HI becomes suppressed by the factor {jj)"^- 

On the other hand, if one is interested in the local 
amplitude of density perturbations for I <^ Aq (which 
is the case if Aq is much greater than the scale of the 
observable part of the universe) one should use Eq. I|12|l . 



IV. CURVATON PRODUCTION DURING 
REHEATING 



The effect which we are going to consider in this sec- 
tion can be very essential, but it is model-dependent. In 
order to make our consideration as simple as possible, 
we will restrict ourselves to the order of magnitude es- 
timates and assume that at the last stages of inflation 
the Hubble constant H, the inflaton mass M and the 
curvaton mass m are roughly of the same order of mag- 
nitude. In reality, we expect that M H ^ m, but 
this does not alter our main conclusion. In fact, it is 
very easy to make our consideration more general and 
precise, following , but we do not want additional 

model-dependent complications to distract us from our 
main qualitative conclusion. 

Our main assumption is that the inflaton field is in 
the hidden sector, and decays to all other fields only due 
to gravitational effects, with the decay rate T (in 
Planck units Mp = 1). We will also assume that the 
curvaton field interacts only very weakly with itself and 
other particles, so that nonrelativistic curvaton particles 
produced as a result of inflaton decay do not thermal- 
ize, and decay at the same time as the classical curvaton 
field and its long-wavelength fluctuations generated on 
inflation. 

At the end of inflation, the energy density of the infla- 
ton field (j) is of order H^, while the energy of the fluc- 
tuations of the curvaton field a is iJ'*. After inflation, 
the energy density of the inflaton is decreasing, while the 
curvaton energy density remains constant until the time 
t ~ m~^, when the long wavelength curvaton field begins 
to oscillate. At that time, the inflaton energy density is 
of order m?. Starting from this moment, the ratio of the 
inflaton and the curvaton energy densities remains con- 
stant, until the inflaton decays. However, when m, Af, 
and H are not different, this is a minor effect. In this 
case, both inflaton and curvaton field begin to oscillate 
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nearly immediately after the end of inflation when their 
energy densities are and H'^ respectively, and their 
ratio, ~ H^^, is preserved until the inflaton decays. 

Ignore for a moment that the inflaton decays to many 
different types of particles and suppose that its energy 
goes entirely to the curvaton particles. Then the energy 
density of the curvaton particles will be greater than the 
energy density of the long-wavelength perturbations of 
the curvaton field by the factor H~^. Note that the con- 
straints on the amplitude of gravitational waves during 
inflation implies that H < 10~^, so that is extremely 
large: H^^ > 10^°, in Planck units. That is why we were 
ignoring numerical factors, the difference between H, M 
and m, etc. The main effect will be present even if we 
relax many of our assumptions. 

Now let us use our estimates for calculating the curva- 
ton contribution to the amplitude of adiabatic perturba- 
tions. This amplitude is given by — , where 5p = rn^aSa. 
As for Per, in all papers on this subject this quantity is 
supposed to be either the energy density of the initial 
classical field a or its generated long wavelength pertur- 
bations. In our scenario the typical value of this last 
contribution is of order (a'^^ ~ H'^. Assuming that it 
dominates and, hence, a ~ /m, we obtain |^ ~ ^, as 
in Eq. \V6l . In order to have ^ ^ 10^^ one needs either 

m ~ \0~^H or, alternatively, we may assume, as many 
people do, that the initial value of the classical curva- 
ton field was much larger than /m, thus, reducing the 
amplitude of perturbations. 

Now let us take into account the curvaton production 
during reheating. First note that the total energy density 
of the produced curvaton particles can exceed the energy 
density of the long-wavelength fluctuations by the factor 
of H^^, and therefore the resulting amplitude of pertur- 
bations is reduced by the factor of H^. Thus even if m is 
comparable with H , we can obtain nevertheless 



(14) 



For H < 10~^, this quantity is smaller than 10^^°. It 
is well below than the level of the standard adiabatic 
perturbations produced by the inflaton field and even 
much smaller than the amplitude of the tensor pertur- 
bations ^ H. In other words, the cr-particles production 
during reheating can dramatically reduce the curvaton 
contribution to the final amplitude of adiabatic density 
perturbations making it negligibly small. 

As we have noted above, one can also decrease the 
contribution of the curvaton perturbations into the final 
amplitude of the resulting adiabatic fluctuations by con- 
sidering the parts of the universe with a very large clas- 
sical curvaton field a. However, if cr > 1, then the field 
a will drive chaotic inflation and will play the role of the 
inflaton instead of the curvaton; see more about it in Sec- 
tion ^1 On the other hand, one can easily check that 



if cr < 1 and one ignores production of curvaton parti- 
cles during reheating, then the amplitude of the curvaton 
contribution — = ~ — is always greater than the 
amplitude of tensor perturbations ^ H produced during 
inflation. 

The possibility of a strong suppression of the curvaton 
contribution due to reheating, Eq. does not mean 

that the final amplitude of the scalar metric perturba- 
tions can be made smaller than the amplitude of the ten- 
sor perturbations. If the contribution of the isocurvature 
curvaton mode is too small, then the resulting amplitude 
of metric perturbations is determined by the usual adia- 
batic perturbations produced by the inflaton field [23 |. 

As we already mentioned, the results obtained above 
are strongly model dependent. One can easily avoid sup- 
pression of the curvaton perturbations by changing the 
assumptions about the mechanism of reheating. The 
main goal of this consideration was to show that in 
the calculations of the curvaton perturbations, in ad- 
dition to the energy of a classical curvaton field and 
its long- wavelength inflationary perturbations one should 
also consider the energy of the curvaton particles pro- 
duced during reheating. This contribution may be ei- 
ther negligible, or extremely important, as in the case 
described above. Since the importance of the effect is 
model-dependent, we will describe it phenomenologically 
by modifying the equation for the curvaton contribution 
to the amplitude of (locally Gaussian) density perturba- 
tions as 



6p_ 

Pa 



2rrra5(j 



aH 



(15) 



Here a is the locally measured amplitude of the classical 
curvaton field and C is a phenomenological parameter 
proportional to the energy density of the curvaton parti- 
cles produced during reheating. 



2^ this equation reproduces 



For C = and CT 
our previous result ^ jj, but as we have seen above 
under certain conditions the constant C may be quite 
large, which leads to a decrease of the amplitude of the 
curvaton contribution to density perturbations, and to a 
different dependence of this amplitude on the local value 
of the field cr. For C ^ 0, the amplitude of the density 
perturbations does not blow up at cr = 0, as before. It 
reaches its maximum at a = C: 



6p 



H 
2^ 



(16) 



One of the interesting consequences of this result is 
that even if one ignores the effects of the curvaton web, 
one can still obtain strongly nongaussian perturbations 
generated by the curvaton mechanism. Indeed, usually 
the amplitude of the locally gaussian fluctuations was 
supposed to be — ~ 10~^. The smallness of this ratio 
implied that the level of nongaussianity was expected to 
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be extremely small (quadratic in ^). However, with an 
account taken of the curvaton particles produced during 
reheating, one can obtain a small amplitude of density 
perturbations even if ^ > 10"^ This may substantially 
increase the level of nongaussianity of adiabatic pertur- 
bations produced by the curvatons. In other words, our 
ignorance of the specific mechanism of curvaton produc- 
tion during reheating makes the degree of nongaussianity 
in the curvaton scenario a free parameter. 

Finally we should note that Eq. (|15(l is valid in a 
more general case, if during the curvaton decay there 
were other particles, giving contribution pp = 2m^C^ to 
the total energy density during the curvaton decay. We 
concentrated on the case where this contribution is dom- 
inated by the curvaton particles, since in this case one 
does not get any undesirable isocurvature perturbations 
if CDM, baryon and lepton asymmetry are generated ei- 
ther by the curvaton decay or after it. One can also avoid 
the isocurvature perturbations even if the energy density 
during the curvaton decay is dominated by radiation, un- 
der the assumption that CDM, baryon and lepton asym- 
metry are produced after the curvaton decay, but not due 
to the curvaton decay . 



V. REHEATING AND THE CURVATON WEB 

The observational consequences of our scenario depend 
on relations between the typical size of the curvaton do- 
main Ao, the size of the observable part of the universe 
Ih and a typical galaxy scale Ig ~ 1Q~^Ih (recalculated 
to the end of inflation). 

As we have seen in the previous section, effects of re- 
heating may allow us to have a small amplitude of lo- 
cally gaussian density perturbations not only for H many 
orders of magnitude greater than m, but also for 
which is only a few times greater than m. This means 
that, depending on the ratio between H and to, a typ- 
ical size of each domain with a positive or negative tr, 

Ao ^ exp^O^^^^ , can take values either much smaller 

or much greater than the scale of horizon Ih- 

We will mention here two interesting regimes: 

If Ih ^ Ig ^ '^Oj the perturbations are strongly non- 
gaussian, which contradicts observational data. On the 
other hand, if the amplitude of these fluctuations is sub- 
dominant compared to the usual adiabatic inflationary 
perturbations, then one can use this regime as providing 
a small controllable amount of nongaussian isocurvature 
and adiabatic perturbations P]. The existence of the 
curvaton web helps to understand the origin of this non- 
gaussianity. 

If Ih ^ Aq, we live in a part of the universe domi- 
nated by nearly gaussian adiabatic perturbations. For a 
proper choice of parameters, such a possibility is consis- 



tent with all available observational data. A nontrivial 
part of this scenario is that the amplitude of adiabatic 
perturbations depends not only on the parameters of the 
model, but also on the local value of the curvaton field 
at the end of inflation, which can vary from to /m 
in different parts of the universe. This helps to justify 
anthropic considerations comparing the properties of the 
universe with different values of the amplitude of density 
perturbations. In this respect, this scenario is similar to 
the scenario studied in [lj| , where the amplitude of den- 
sity perturbations could take different values in different 
parts of the universe due to inflationary perturbations of 
the Brans-Dicke field. Recently this possibility was stud- 
ied in in relation to the anthropic principle and string 
landscape scenario. We will return to the discussion of 
this issue in Section IVIII The picture of the curvaton 
web allows not only to discuss various probability distri- 
butions, but also to have a visual understanding of the 
structure of the universe in this scenario. 

The details of this picture depend on creation of a- 
particles during reheating discussed in the previous sec- 
tion. In the case C = considered in Section IIIII (no 
curvaton particle production during reheating), the max- 
imum of density perturbations |^ is reached near the 
walls cr = separating the exponentially large domains 
with positive and negative values of a. However, for 
C 7^ 0, the maximum would correspond to a finite value 
\<Jm\ = C (see Eq. ^^). If C > H^/m, then the ampli- 
tude of density perturbations will take its largest values 
deep inside these domains, when a ^ /m, and it will 
be minimal near the walls with a ~ Q. On the other 
hand, if C < /m, then the amplitude of density per- 
turbations will be small inside the domains, it will grow 
on the way towards the walls with a ~ Q, but then it will 
reach some maximal value and vanish near the walls with 
(7 = 0. This means that the places where the amplitude 
of density perturbations is maximal form two walls sur- 
rounding from each side the wall with cr = 0. Thus the 
effects of reheating may affect in a rather nontrivial way 
the structure of the curvaton web. 

Our results are illustrated by a set of two figures, which 
show the spatial distribution of the square of the ampli- 
tude of the curvaton density perturbations based on com- 
puter simulations of the curvaton web for C < H^/m, see 
Fig. 121 and Fig. 13 These two figures show why did we call 
the resulting structure 'the curvaton web.' 

These figures clearly explain the main source of non- 
gaussianity in the curvaton scenario of Ref. ^2,]. They 
also suggest a rather exotic possibility which may occur 
if we live not very far from the walls shown in Fig. |21 
In this case the amplitude of density perturbations will 
be slightly different at the two opposite parts of the sky, 
depending on the distance from the nearby wall. 
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FIG. 2: The spatial distribution of the square of the amplitude 
of density perturbations near the walls (7 = 0, for C < H^/m. 
The amplitude of density perturbations vanishes at the walls 
a = 0. These walls are surrounded from both sides by the 
walls where the amplitude of density perturbations reaches 
its maximal value. 




FIG. 3: Contour plot for the square of the amplitude of the 
density perturbations shown at the previous figure. White 
regions correspond to vanishing amplitude of density pertur- 
bations, black regions correspond to maximal absolute val- 
ues of density perturbations. This figure confirms that the 
two walls with large amplitude of density perturbations (see 
Fig-HJ surround the domain boundaries with a — where the 
amplitude of density perturbations vanishes. (At the previous 
figure, the regions with a — were hidden by the walls.) 



VI. CURVATONS AND ETERNAL CHAOTIC 
INFLATION 



Curvatons and chaotic inflation 



Until now we were assuming that the Hubble constant 
and the mass of the curvaton field remain constant dur- 
ing inflation. This possibility can be realized, for ex- 
ample, in new or hybrid inflation. In chaotic inflation 
the Hubble constant changes signiflcantly during infla- 
tion. Moreover, the curvaton mass can depend on the 
curvature, = tuq + aH^, see 0. In this Section we 
consider how the results obtained in the previous sec- 
tions are modified in this case. We will also assume, as 
in the previous sections, that initially the scalar field a 
was small; for a discussion of a more general situation see 
[Tol |25| and also the end of section IVIIjl 

We consider the model where inflation is driven by 
the field with an effective potential F(0) and assume 
that the curvaton mass is much smaller than the Hubble 
constant during inflation, that is <C H^- In this case 
the long-wavelength modes of the curvaton field are in 
the slow-roll regime and satisfy the equation 



3H&k + m ak — 



(17) 



and, hence, for every mode with a given comoving 
wavenumber k we have 



dt m^^- ^^^^ 

Integrating over k and taking into account the contri- 
bution of the new ly g enerated from quantum fluctuation 
modes we obtain 



2m^ 



47r2 



(19) 



dt 3H 

The solution of this equation with the initial condition 



is 



1 

4^ 



(t) exp 



ti 



Hit) 



dt\ di . (20) 



If H and m are both constant we obtain from here the 
familiar result as i ^ oo. In the case of chaotic in- 
flation it is convenient to rewrite this equation using the 
value of the inflaton fleld as a time variable instead of 
t. Taking into account that 

H^^^, (21) 

dt m' 3 ' ^ ' 

where we again use the system of units Mp — 1 , equation 
(Pn)l becomes 



(0)) 



1 

12^ 



exp 



m 



(22) 
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where (t)i is the initial value of the inflaton field. Let us 
consider separately two cases: vn? — const and = 
aH^, where a is some small constant, a < 1. 



Curvaton with m = const 



We begin with simple inflation due to the massive in- 
flaton field with potential V — iM^0^. In this case Eq. 
(|22|l implies that at the end of inflation, when ~ 1, we 
have 



(4 



2m^ 
A/2 



In particular, for ^ A^P one obtains [2^ 



M 



2a4 



1927r2 



(23) 



(24) 



These fluctuations lead to the same consequences as a 
classical scalar field c which is homogeneous on the scale 
and which has a typical amplitude 



^0 = VF> 



(25) 



As we see, the value of <tq depends on the initial value 
of the field </>. This result has the following interpreta- 
tion. One may consider an inflationary domain of initial 
size H~^(<f>i). This domain after inflation becomes expo- 
nentially large. In the model with V{(1>) = -^0^ its size 
becomes 



largi 

fir 



i)exp 



(26) 



One can easily check that the largest contribution to (o-^) 
is given by perturbations with the wavelength 0(Ao). 



If one takes 0; ^ 0(1/ v M), a typical amplitude of 
the field a becomes much greater than 0(1) in Planck 
units, which means that this field becomes capable of 
supporting chaotic inflation.^ One may call this effect 
'the curvaton- inflaton transmutation.' 

In this regime the field a typically serves as an infla- 
ton rather than as a curvaton. Indeed, in this case we 
consider two fields, ip and a, in a universe with a Hubble 
constant determined by a total energy density of both 
fields. Their combined dynamics can be easily under- 
stood in the approximation H const. In this case 



Ob- 



4> — (/'icxp (^—^jj^^, whereas a — do exp 

viously, fox m <^ M < H the field a practically does 
not change until inflationary stage driven by the field (j) 
is over. After that, the field a drives subsequent stages 
of infiation, producing usual infiaton perturbations with 
amplitude 0{m). 

Thus, if the field (f> was originally greater than 
0{1/\/M), then the last stage of infiation is driven by 
the lightest field tr. This happens even if the initial value 
of a is smaller than the Planck value; see |23] for a related 
discussion. 

Until now we were assuming that the scalar field (f> 
obeys classical evolution described by Eq. H21|l . However, 
if <j)i > 0{l/y/M) then inflation is eternal [HE^. Thus 
in eternal inflation scenario based on the model with two 
fields, the last stage of inflation is typically driven by the 
lightest field. In other words, the lightest field a typically 
is not a curvaton but an inflaton. 

Moreover, even if eternal inflation is initially supported 
only by the field 0, in some parts of the universe the 
quantum fluctuations bring this fleld higher and higher, 
up to the Planck density. In those regions, the Hubble 
constant grows up to the Planck value H ^ 1, and (o-^) 
eventually becomes 

{a') ^ Oim-') . (27) 



The energy density of the fluctuations of the field a also 
approaches Planck density: im^(tT^) = 0(1). This is 
not surprising; in fact this always happens in the eternal 
inflation scenario based on the simplest chaotic inflation 
which allows the nearly Planckian energy density 'T8','28j. 
But this means that as a result of quantum fluctuations 
produced during eternal inflation, the scalar field a typ- 
ically becomes much greater than the field (j). Indeed, 
the field a can take any value in the interval \a\ < 1/m, 
which is much greater than |(/)| < 1/M, in Planck units. 
Later on, the field (j> rolls down faster, and the last stage 
of inflation is eventually driven by the lighter field a. 

How these results depend on the choice of the inflaton 
potential V (</))? Assuming that during inflation driven 
by the fleld (j) the curvaton mass m is negligibly small 

compared to the Hubble constant, and the ratio ^-^^f'-* 
grows at large (j>, we obtain from H22|) the following esti- 
mate for (a^) at the end of the inflation: 



V (0.) 



(28) 



Therefore if at the beginning of inflation (i.e. at = (^i) 
the condition 



^ A closely related thought somewhat different possibility was 
mentioned by David Wands in his talk at the workshop after 
the conference "The Next Chapter in Einstein's Legacy" at the 
Yukawa Institute, Kyoto. 



(</)) </) > y (0) 



(29) 



is satisfied, then (u^) > 1, and the curvaton a becomes 
an infiaton. It is easy to verify that, for example, for 
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the power-law potential V 
reproduction, 



Xcj)", the condition of self- For a ^ (j) ^ we have 



(30) 



coincides with (|29|l . Therefore all results which we de- 
rived for the massive inflaton potential are valid in this 
case as well. 

Thus we are coming to the conclusion that if eternal 
inflation takes place in our model (or if the field a was 
greater than 0(1) from the very beginning), then the 
scalar field cr typically becomes extremely large and plays 
the role of the inflaton rather than the curvaton. We will 
return to this conclusion and its consequences in Section 

Ivlil 



Curvaton with m = aH , a < 1 



Another interesting case is when the curvaton mass 
squared during inflation is proportional to the curvature 
R — \2H^ Mass terms ~ aH^ appear in many 

models of supergravity, and also in the models of the 
field a non-minimally coupled to gravity. Substituting 
w? = aH"^ = f in (123), we obtain 



1 



127r2 



(*) 



exp 



/ 



V 



2a f V 



3 J V 





-d(j) d(f) 



(31) 

For the case of the power-law potential V the integral in 
(pTT|l can be calculated exactly. In particular, for V = 
iM^0^, we obtain 



loTT^a \ a \ a J 



In the case a (p, ^, this expression simplifies to 



192^ 



(32) 



(33) 



By taking (f) <^ (pi, we reproduce our previous result 124(1 
for the amplitude of perturbations of the field tr at the 
end of infiation driven by the field (f). Therefore, for 
a ^ M, the mass term aH^ is irrelevant at the self- 
reproduction scale (pi ~ M~^/^ and below. As a result, 
in this regime the curvaton a typically becomes an infla- 
ton and drives the last stages of inflation. 

On the other hand, for sufficiently large a, the behavior 
of the curvaton filed in chaotic inflation is well described 
by the model with constant m and H considered in Sec- 
tion ^1 Indeed, if a ^ ^7^' then, as soon the held <f) 
becomes much smaller than <j)i, one can neglect the last 
term in 132II . and the result becomes 



M 



/f2 



IGTT^a 



(34) 



M2,/)2 3ij4 (0) 



IGtt^q; Sn'^m? {(j)) 



(35) 



Thus, the distribution of the flcld a behaves as a Bunch- 
Davies distribution for a field with an adiabatically de- 
creasing mass m? k, aH^. A typical dispersion of the 

field cr slowly decreases, remaining of the order ^^—^ ~ 

This result is easy to understand. As it was shown 
in Sectional the Bunch-Davies distribution |j3Jl is estab- 
lished during the time t ^ All results obtained in 
Section^will remain valid at all stages of inflation where 
the Hubble constant did not change much during the time 
i ~ ^ = One can show that this is equivalent to 
the condition 



a > 



(36) 



This adiabatic regime ends when 0^ becomes smaller 
than a~^. Thus, at the last stages of inflation Eq. 
yields 



3M2 



(37) 



Note that if a ^ M, then (cr^) < 1, and the curvaton 
does not become an inflaton. 

As we have shown in Q, the slope of the spectrum of 
curvaton perturbations produced in this model depends 
on the parameter a. If a <t; 1, the spectrum will be nearly 
flat, and perturbations will be locally gaussian. If a is 
sufficiently large, then the typical size of the cells of the 
curvaton web become small, the perturbations become 
strongly nongaussian, and they may have blue spectrum. 

We should emphasize again that the curvaton pertur- 
bations appear in addition to the usual adiabatic pertur- 
bations produced by the inflaton fluctuations. By com- 
bining these two contributions, one can obtain adiabatic 
perturbations with a controllable degree of nongaussian- 
ity. Also, since each of these two contributions has its 
own spectral index, a combination of these two spectra 
may look like a spectrum with a running spectral index. 



VII. CURVATON WEB AND STRINGY 
LANDSCAPE 



In this section we will discuss an additional reason to 
be interested in the curvaton models and the curvaton 
web. 

With the development of inflationary scenario it be- 
came clear that the universe may consist of exponentially 
large number of exponentially large locally homogeneous 
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domains of different types, which allowed us for the first 
time to obtain a scientific justification of anthropic prin- 
ciple m m m lal E3- T^s possibility became espe- 
cially interesting when it was realized that string theory 
may have an enormously large number of different vacua 

One may wonder where is our place in the vast stringy 
landscape? The first thing to say is that we simply can- 
not live in most of these vacua. For example, we cannot 
live in a non-compactified lOD space even if its total vol- 
ume is extremely large. This is a simple superselection 
rule. One may try to go further and evaluate the proba- 
bility of emergence of life in a world with different masses 
and coupling constants. This is what was traditionally 
done by those who explored various consequences of an- 
thropic principle prior to invention of inflationary cos- 
mology. 

Inflation allows us to look at this problem from a dif- 
ferent perspective. One may assume that we are typical 
observers, and the total number of observers is propor- 
tional to the total volume of the universe. Therefore the 
probability to live in a part of the universe with given 
properties depends not only on these properties, but also 
on the volume of this part of the universe as compared to 
the volume of other parts. Since volume is exponentially 
sensitive to parameters of inflationary theory and ini- 
tial conditions for inflation, investigation of the volume- 
weighted probability distribution of different parts of the 
universe can give us a very powerful tool for understand- 
ing of our place in the universe. 

The volume-weighted probability distribution in infla- 
tionary cosmology was introduced in ,28^ ,36J . A detailed 
investigation of this probability distribution and its ap- 
plications to anthropic principle was performed in 'l8l, 
which was followed by many other works on this subject, 
see e.g. [sOi 0, 113, S| and references therein. 

There are many problems related to the use of this 
probability distribution, including the choice of a proper 
probability measure and a general issue of its interpreta- 
tion. These issues could seem abstract and academical, 
but recently the situation changed. After the invention 
of the mechanism of moduli stabilization in string theory 
and the development of the concept of string theory land- 
scape ^2, J3., 34^ 35} the problem of the choice of a most 
probable vacuum state became one of the most important 
problems of modern physics. Many string theorists and 
cosmologists began actively debating this question. It is 
not our goal to discuss all of these issues in this paper. 
Here we will only mention one aspect of this problem 
related to the curvaton web. 

To explain the basic issue, let us consider, following 
p^. I30II39I I . the simplest class of inflationary models with 
a potential = ^^(f^, which we discussed in Sect. 

IVII During the rolling of this field from its initial value (f>i 
the universe expands as e'^^^^, see Eq. H26|) . Suppose in- 



flation began at the Planck boundary V((/)i) — ^(/>^ — 1. 
Then, ignoring the issues of eternal inflation, the volume 
of the universe during inflation in this scenario grows 
exp ( 21^^) times. If instead of that we assume that infla- 
tion begins at the eternal inflation boundary ^ 1/ ^/M, 
one finds that the volume of the universe grows by a fac- 
tor of exp (0(M~^)). In either case, the total volume of 
the universe after inflation grows exponentially when AI 
decreases. 

Thus if the inflaton mass M can take different values 
in different parts of the universe, one may argue that 
the parts with a smaller M will be exponentially bigger, 
so they will contain exponentially greater number of ob- 
servers, and therefore a typical observer will live in a part 
of the universe with the smallest possible inflaton mass. 

On the other hand, the amplitude of density perturba- 
tions produced in this scenario is 0, 0, 0| 



Sp 
P 



y3/2 
V 



M(j3^ - M 



(38) 



at the end of chaotic inflation, where = 0(1). Since 
the volume of the parts of the universe produced by in- 
flation is exponentially sensitive to the value of M , one 
may argue that it would be exponentially more proba- 
ble to live in the universe with the amplitude of density 
perturbations much smaller than its present value. 

This is a serious problem, which was especially clearly 
described in Ref. |2J| on the basis of a combined anal- 
ysis of the probability distribution for the cosmological 
constant and the amplitude of density perturbations. 

One should note that this problem appears mostly be- 
cause of the assumption that the probability to live in a 
given part of the universe is proportional to the degree 
of inflation. While this is a reasonable assumption [l8l| . 
its validity is not absolutely clear; see e.g. a discussion 
in [3lll40||. In particular, if one uses the probability mea- 
sure associated with the volume at the hypersurface of 
time measured in units of , which was one of the two 
possible measures advocated in 18], then the resulting 
probability distribution does not contain any factors of 
the type of e*^'/^. In addition, we do not actually know 
whether it makes any sense to compare different causally 
disconnected parts of the universe with different param- 
eters |3l) |. 

These issues deserve a separate investigation, which 
goes beyond the scope of the present paper. Still, with 
all of these uncertainties in mind, the problem described 
above should not be discarded. In order to resolve this 
problem, it was suggested in [s^l that the cosmological 
perturbations of metric should be produced not by the 
inflaton field but by cosmic strings, in which case the 
amplitude of perturbations does not depend on M and 
on the duration of inflation. This possibility by now is 
ruled out by cosmological observations. Therefore in |23| 
it was proposed to use the curvaton mechanism of gen- 
eration of density perturbations. The main idea is that 
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in the simplest versions of this scenario ^ see 
Eq. (P). This amphtude does not depend on the mass 
of the inflaton field, so one can maximize the amount of 



inflation by minimizing M , without affecting 



Sp 



One of the problems with this proposal is that in this 
scenario the probability to live in the universe with the 
amplitude of density perturbations greater than some 
given value and with the anthropic constraints im- 
posed on the combination of ^ and the cosmological 

constant A, is proportional to 24] . This result sug- 

gests that we should live in the universe with ^ 3> 10 "'^j 
which also disagrees with the observational data. In order 
to avoid this problem it was suggested to consider more 
complicated, multi-component curvaton models . 

In this respect, our investigation of the curvaton web 
brings both good and bad news. The good news is that 
if one takes into account production of the curvaton par- 
ticles during reheating, then the amplitude of the cur- 
vaton perturbations does not blow up at small a. In- 



stead it has a maximal value given by 



Sp 



H 
27rC' 



see 

Eq. H16|l . This would be consistent with observations for 
C - IQ^H. 

The bad news is that this idea requires a specific care- 
fully tuned mechanism of reheating after inflation. More- 
over, this idea does not work at all in the simplest cur- 
vaton models. For example, as pointed out in Section 
IVll in the eternal inflation scenario in the model with 
the curvaton potential m?a'^ /2 and the inflaton poten- 
tial M'^cjP /2, the field a typically plays the role of the 
inflaton instead of the curvaton. 

One can avert this conclusion by considering the field 
a with the mass m? — rriQ + aH^ with a sufficiently 
large a, as we did in [9|; see Section IVl CI But this 
makes the model more complicated and prompts to look 
for alternatives ideas. 



The simplest way to obtain 



In + aH'^ is to add 



the term — ^ct i? to the Lagrangian. But if one does 
it for the curvaton, one may do it for the inflaton as 
well. By adding the term — to the Einstein La- 
grangian ■^R, one would make it impossible for the field 
(j> to be greater than 1 / because the effective gravita- 
tional constant G — (87r(l — ^0^))^^ becomes singular at 
(j) = l/\/J- In this model the maximal amount of infla- 



tion can be estimated by e' 



4't/i 



= 1/4? 



This quantity no 



longer depends on A/, so the whole problem discussed in 
^^.3^,32,] disappears. 

More generally, thisproblem arises only because the 
models considered in 0, 113, H^l depend on a single pa- 
rameter, such as the inflaton mass M, or a coupling con- 
stant A in the potential As we just mentioned, this 
problem disappears if we add a second parameter, ^, and 
keep it fixed. Of course, one may argue that this param- 
eter must be minimized in its turn, in order to maximize 



the degree of inflation. But with each such step the reli- 
ability of the basic arguments becomes less certain. 

In this respect, it is quite instructive to consider a 
different modification of the theory. Let us add a con- 
stant Vq to the potential M^^^/2, and assume that in- 
flation ends when the field reaches some point 0o with 
M^(/)o/2 <^ Vb, as it happens in the hybrid inflation sce- 
nario 1^. In this case, as before, the total duration of 
inflation will be determined by the behavior of the uni- 
verse at extremely large 0, where M'^(jP/2 '--^ 1 ^ Vq. 
This duration will be almost the same as in the simple 
model M'^cjP /2. Meanwhile the amplitude of density per- 
turbations during the last 60 e-folds will be determined 
by inflation in the re gim e 

< Vo, and it will be 

proportional to A/~^ uM- 



Sp 
P 



v 



3/2 



(39) 



In this case one would have another problem: Long dura- 
tion of inflation, as before, favors small M, but now small 
M lead to unacccptably large density perturbations. 

So what is the problem? Is it an unacccptably small 
density perturbations, as in the model M^^^S, 

or an 

unacccptably large density perturbations, as in the model 

M2(/)2/2 -I- Vo? 

The root of this ambiguity is easy to understand. The 
degree of inflation in this scenario is determined by the 
shape of the potential at extremely large (f>, whereas the 
amplitude of density perturbations is determined by the 
shape of the potential at relatively small (j). The problem 
discussed in [2J, ISfl |33 appears only if one believes that 
the shape of the inflaton potential is determined by a sin- 
gle parameter and does not change during an incredibly 
large interval of variation of the inflaton field. We do 
not think that this requirement is generically satisfied by 
inflationary models in stringy landscape. 

The simplest way to address this problem is to max- 
imize the total amount of inflation by properly tuning 
the part of the potential describing the main part of the 
process of inflation (which corresponds to (j) ^ 0(10) in 
the simplest chaotic inflation models) , and then tune the 
part describing the last 60 e-folds of inflation by tuning 
the remaining part of the potential {4> ^ 0(10) in our ex- 
ample). As our example shows, one can easily tune the 
amplitude of perturbations of metric anywhere in the in- 
terval from to 0(1) without altering the early stage 
of inflation which is responsible for the total degree of 
inflation. 

Therefore we believe that it would be somewhat pre- 
mature to use anthropic considerations for justification 
of the curvaton scenario. 
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VIII. CONCLUSIONS 

In this paper we discussed the global structure of the 
universe in the simplest curvaton scenario proposed in 
our paper Q. This scenario is based on the theory of 
two massive scalar fields, the inflaton and the curva- 
ton cr, interacting with each other only gravitationally. 
We have shown that stochastic processes during eternal 
inflation in this scenario divide the universe into many 
exponentially large domains containing either positive or 
negative values of the curvaton field a. The amplitude of 
the curvaton perturbations in this scenario depends on 
the process of reheating after inflation and has a rather 
nontrivial behavior as a function of the distance from 
the walls a ~ separating the domains with positive 
and negative a. We called the resulting structure, which 
was shown in Figs. 13131 'the curvaton web.' 

Depending on the parameters of the model and on our 
own position with respect to the curvaton web, we may 
live in a universe with strongly nongaussian adiabatic 
perturbations, or in the universe with locally gaussian 
adiabatic perturbations with the amplitude depending on 
our position, or in the universe dominated by the usual 
gaussian inflaton perturbations with a small admixture 
of nongaussian curvaton perturbations. In the last case, 
one may relax many constraints imposed by the require- 
ment of a secondary reheating and creation of dark mat- 
ter and baryons after the curvaton decay, and instead 
of the adiabatic perturbations produced by the decaying 
curvaton study a small admixture of isocurvature non- 
gaussian perturbations. Moreover, in different parts of 
the universe the same scalar fields may play either the 



role of the curvaton or the role of the inflation; the effect 
which we called 'the curvaton-inflaton transmutation.' 

This example shows that in addition to a discrete spec- 
trum of possibilities, which appear after compactification 
in string landscape scenario, one may have a continuous 
spectrum of possibilities, which appear because of a com- 
plicated dynamics in the early universe. The availabil- 
ity of so many new possibilities is encouraging, but also 
rather disturbing. The curvaton scenario is much more 
complicated than the usual scenario where the adiabatic 
perturbations are produced by the inflaton fluctuations, 
so we believe that it should be used sparingly. On the 
other hand, in certain situations the curvaton scenario 
can be quite useful. In particular, it provides by far the 
simplest way to obtain a controllable amount of nongaus- 
sian perturbations of metric 9] . It has many other inter- 
esting and unusual features, some of which were discussed 
in this paper, so it certainly deserves further investiga- 
tion. 
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